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Abstract 

The dynamics of a complex system is usually recorded in the form of time series, which can be studied through its 
visibility graph from a complex network perspective. We investigate the visibility graphs extracted from fractional 
Brownian motions and multifractal random walks, and find that the degree distributions exhibit power-law behaviors, 
in which the power-law exponent a is a linear function of the Hurst index H of the time series. We also find that the 
degree distribution of the visibility graph is mainly determined by the temporal correlation of the original time series 
with minor influence from the possible multifractal nature. As an example, we study the visibility graphs constructed 
from two Chinese stock market indexes and unveil that the degree distributions have power-law tails, where the tail 
exponents of the visibility graphs and the Hurst indexes of the indexes are close to the a ~ H linear relationship. 

Key words: Visibility graph; complex networks; power-law distribution; fractional Brownian motion; multifractal 
random walk 

PACS: 89.75.Hc,05.40.-a, 05.45.Df, 05.45.Tp 



1. Introduction 

Complex systems are ubiquitous in natural and social sciences, where the constituents interact with one another 
and form a complex network. In recent years, complex network theory has stimulated explosive interests in the 
study of social, informational, technological, and biological systems, resulting in a deeper understanding of complex 
systems d S S 01 ■ However, for many complex systems, it is hard to obtain detailed information of interacting 
constituents and their ties, which makes the underlying network invisible. Instead, we are able to observe and record a 
time series generated by the system. For such cases, time series analysis becomes a crucial way to unveil the dynamics 
of complex systems. There are also some efforts to map time series into graphs to study time series from the network 
perspective, which amounts to investigating the dynamics from the associated network topology. 

For a pseudoperiodic time series, one can partition it into disjoint cycles according to the local minima or maxima, 
and each cycle is considered a basic node of a network, in which two nodes are deemed connected if the phase^space 
distance or the correlation coefficient between the corresponding cycles is less than a predetermined threshold fl3|] . We 
note that a weighted network can also be constructed if the phase space distance or the correlation coefficient is treated 
as the weight of a link. This method for pseudoperiodic time series can also be generalized to other time series, where 
a node is defined by a sub-series of a fixed length as the counterpart of a cycle, which has been applied to stock prices 
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Another method for network construction from time series is based on the fluctuation patterns Htl- In this 
approach, each data point is encoded as a symbol R, r, D, or d, corresponding to big rise, small rise, big drop 
and small drop, respectively. The time series is then transformed into a symbol sequence. Defining a fluctuation 
pattern as an n-tuple consist of a string of n symbols, the symbol sequence can be further mapped into a sequence 
of non-overlapping n-tuples. The n-tuples are treated as nodes of the constructed network. Therefore, the number 
of nodes does not exceed n A . Two nodes are connected if the associated n-tuples appear one after the other in the 
pattern sequence. Furthermore, the edge weight between two nodes can be defined as the occurrence number of two 
successive patterns in the sequence. This approach has been applied to study the price trajectory of Hang Seng index 
01. 

A third method is to convert time series into visibility graphs based on the visibility of nodes JjJ. For simplicity, 
consider an evenly sampled time series {y t : t = 1,2, ■ ■ ■ ,N}. Each data point of the time series is encoded into a node 
of the visibility graph. Two arbitrary data points y, and yj have visibility if any other data point located between 
them fulfills 

yj-yk ; yj-yi 
j-k j-i 

Two visible nodes become connected in the associated graph. An example of a time series containing 16 data points 
and the associated visibility graph derived from the visibility algorithm is illustrated in Fig. [T] By definition, any 
visibility graph extracted from a time series is always connected since each node sees at least its nearest neighbor(s) 
and the degree of any node y, with 1 < t < N is not less than 2. In addition, a periodic time series converts into 
a regular graphs, whose degree distribution is formed by a finite number of peaks related to the series period, while 
random time series lead to irregular random graphs @]. It is also found that visibility graph is invariant under affine 
transformations of the series data since the visibility criterion is invariant under rescaling of both horizontal and 
vertical axes, and under horizontal and vertical translations 0. 




Figure 1 : Example of a time series containing 1 6 data points (upper panel) and the associated visibility graph derived from the visibility algorithm 
(lower panel). 



Degree distribution p(k) is one of the most important characteristic properties of complex networks H10Q . The 
degree distribution of the visibility graphs of several specific examples of time series have been investigated |3]. For 
a random time series extracted from an uniform distribution in [0,1], the degree distribution of the visibility graph 
has an exponential tail p(k) ~ e~ k ^ . Alternatively, the visibility graphs of Brownian motions and Conway series are 
scale-free, characterized by a power-law tail in the degree distribution: 



p(k) ~ k~ 



(2) 



where a = 2.00 + 0.01 for Brownian motions and a = 1.2 + 0.1 for Conway series. It is also conjectured that 
the temporal correlation of the time series (characterized by its Hurst index H) might have influence on the degree 
distribution of its visibility graph |9]. In this work, we test this projection based on extensive numerical simulations. 
Specifically, fractional Brownian motions (FBMs) ll 111 and multifractal random walks (MRWs) fl2| are synthesized 
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to investigate the influence of autocorrelation and multifractality on the degree distribution. 



2. Numerical analysis 



2.1. Generating FBMs and MRWs 

There are many different algorithms for the generation of fractional Brownian motions 1 1311 and we adopt a 
wavelet-based algorithm to simulate FBMs [14]. On the other hand, a multifractal random walk can be generated 
by the cumulative summation of the increments 



Ay, 



(3) 



where e, is a fractional Gaussian noise with Hurst index H{ n , to, is a correlated Gaussian noise, and they are independent 
1 1211 . We use the detrended fluctuation analysis IU51I16I1 to verify if the resultant Hurst index of the generated signals is 
identical to the input value of H; n in the algorithms. For each H m , we generate 10 realizations and calculate the mean 
Hurst index H. The results are presented in Fig. [2] The top panel of Fig.|2]shows that the estimated Hurst indexes of 
the synthesized FBMs are very close to the input value H m with minor deviation for small H m . For MRWs, we find 
that H = H m when H\ n ^ 0.5 and a systematic deviation H > H m when Hi n ^ 0.5. In addition, we have confirmed 
that the generated MRW signals possess multifractal nature. 
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Figure 2: (color online.) Dependence of the Hurst index H of the simulated (a) FBMs and (b) MRWs, determined by detrended fluctuation analysis 
on the input Hurst index H m in the two synthesis algorithms. 



2.2. Numerical results 

We have investigated FBMs with the input Hurst index H m ranging from 0. 1 to 0.9 in the spacing of 0. 1 . For each 
H m , we repeat the simulation 100 times and each simulation gives a FBM signal with the size N = 50, 000. For each 
FBM signal, a visibility graph is constructed and its empirical degree distribution is determined. We find that the 100 
distributions almost collapse onto a single curve. This enables we to put all the data of the 100 graphs to construct 
the empirical degree distribution to gain better statistics. Three typical empirical degree distributions of the visibility 
graphs converted from FBM series with different Hurst indexes are depicted in Fig.[3ja). Nice power-law behaviors 
are observed in the distributions, followed by faster relaxation. We note that the visibility graphs of other FBMs also 
exhibit power-law tails in the degree distribution. 

The situation is very similar for the MRW case. We have simulated MRW signals of size N = 50, 000 with 
different input Hurst index H m ranging from 0.05 to 0.95 with an increment of 0.05. For each H m , 100 MRW signals 
are simulated and then converted to 100 visibility graphs. Three typical empirical degree distributions of the visibility 
graphs are depicted in Fig.[3jb) for different Hurst indexes H (not Hm). All the distributions exhibits nice power laws 
with faster decay for large degrees. 
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Figure 3: (color online.) (a) Empirical degree distributions of the visibility graphs converted from MRW series with different Hurst indexes 
H = 0.33, 0.5 and 0.8. (b) Empirical degree distributions of the visibility graphs converted from FBM series with different Hurst indexes H = 0.2, 
0.5 and 0.8. 



The power-law exponents a of the distributions are calculated in the scaling ranges. Figure|4]shows the dependence 
of the power-law exponents a on the Hurst indexes H. Both curves show a nice linear relationship: 

a(H) =a-bH . (4) 

A least-squares regression gives a = 3.35 and b = 2.87 for FBMs and a = 3.19 and b = 2.55 for MRWs. We find that 
the multifractal nature of the MRWs has minor influence on the degree distributions of the visibility graphs. 
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Figure 4: (color online.) Dependence of the power-law exponent a on the Hurst index H for (a) FBM visibility graphs and (b) MRW visibility 
graphs. The straight lines are the least-squares fits of Eq. (3). 



We note that, the linear relationship between the tail exponent and the Hurst index was also found for fractional 
Brownian motions independently Jl7ll . 



3. Applications to financial data 

In this section, we apply the visibility graph method to financial data. Specifically, two indices of the Chinese 
stock market are considered. The organized stock market in mainland China is composed of two stock exchanges, the 
Shanghai Stock Exchange (SHZE) and the Shenzhen Stock Exchange (SZSE). Shanghai Stock Exchange Composite 
index (SHCI) and Shenzhen Stock Exchange Component index (SZCI) are the representative indices for the two stock 
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exchange, respectively. Our analysis is based on the 1-min data of the two indices. The time series spans from 2 
January 2001 to 28 December 2007 for the SHCI and from 4 January 2002 to 28 December 2007 for the SZCI. The 
temporal evolution of the price trajectories of the two indices are illustrated in Fig. [5] According to Fig. [5] the Chinese 
stock market was in a bearish phase from 2001 to 2005, which is known as an antibubble 11811 . and then the market 
reversed and produced a very marked bubble which burst at the end of 2007. In order to test if the market phase has 
impact on the results, we partition each time series into two sub-series delimited on 3 1 December 2005, corresponding 
to the bear period and the bull period, respectively. 
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Figure 5: Temporal evolution of the price trajectory of the Shanghai Stock Exchange Composite index and the Shenzhen Stock Exchange Compo- 
nent index. 



For each index, we obtain three visibility graphs corresponding to the bear period, the bull period, and the whole 
time period. The empirical degree distributions of the three visibility graphs are determined. Fig.|6]shows the results 
of SHCI. We note that the results for the SZCI are very similar. 




Figure 6: (color online.) Empirical degree distributions of the three visibility graphs of the Shanghai Stock Exchange Composite index. 



We observe that all the distributions have heavy tails and the degree distribution in the bear period can be well 
modeled by a power law (O in the right part. The rapid decay does not mean that the power law is truncated but rather 
reflects the fluctuations at finite size II 1 911 . The tail exponent is estimated to be or = 1.95. For the case of the bull 
period, there is a hump in the tail. This phenomenon is more evident for the whole period, which is caused by the fact 
that there are more points in the bull period that can see more previous points (large degrees). There are two power 
laws in the degree distributions for the whole time period case, illustrated by two parallel dashed lines in Fig. [6] The 
estimates of the tail exponents a for SHCI and SZCI are listed in Table [T] 

In order to compare the empirical results with the numerical results in Section|2] we determine the Hurst indexes 
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Table 1 : Comparison of the estimated tail exponents a and the "predicted" tail exponents a' . The last row is e = (a - a' )/a' 
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1.67 
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a' 


1.85 


1.92 
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1.88 
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0.05 


-0.02 


0.06 


0.05 


-0.11 


-0.07 



of the log returns for each index by applying the detrended fluctuation analysis. The results are listed in TableQ] One 
finds that the Hurst indexes are undistinguishable from H = 0.5. The multifractal detrended fluctuation analysis I20I1 
confirms that all the time series possess multifractal nature. Therefore, the tail exponents can be "predicted" according 
to the linear relationship a' = 3.19-2.55//. The predicted tail exponents are also presented in TableQ] It is found that 
the discrepancy between a and a' is not large, which is quantified by the relative difference e - (a - a') I a' shown in 
TableQ] Fig. Q] further illustrate this point by the scatter plot of the six data points (//, a). These data points are very 
close to the dot-dashed line a' — 3.19 - 2.55//, showing that the empirical results are consistent with the numerical 
analysis in Section|2] 
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Figure 7: (color online.) Scatter plot of (H,a) of the six time series comparing the empirical results with the numerical results predicted by 
a' = 3.19 - 2.55// (dot-dashed line). 



4. Conclusion 

In summary, we have studied the degree distributions of visibility graphs extracted from fractional Brownian 
motions and multifractal random walks. We found that the degree distributions exhibit power-law behaviors, in which 
the power-law exponent is a linear function of the Hurst index inherited in the time series. In addition, the degree 
distribution of the visibility graph is mainly determined by the temporal correlation of the corresponding time series, 
and contains minor information about the multifractal nature of the time series. The linear relation (0| provide a 
possible tools for the determination of H of a time series from its visibility graph However, cations should be 
taken since the increments distribution of the time series might also have impact on the degree distribution. 
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